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' ■ Abstract 

^ ' In this paper we present an exact solution of Einstein's field equations, de- 

^ . scribing a rotating stationary axisymmetric dark matter, which is an extension of 

' non-rotating dark matter solution of Paper I. We find that the rotating stationary 

solution is a non- asymptotic, Pctrov type D space-time whose energy-momentum 
(N ■ tensor admits a dark fluid having negative pressure and the energy equation of 

state parameter with minus sign. We also find that the time-like vector fields 
of the matter distribution is expanding, shearing and rotating with acceleration. 
. It is also found that, due to the negative pressure, the energy-momentum tensor 

T-H ' of the solution violates the strong energy condition leading to the repulsion of 

. the gravitational field of the rotating space-time geometry. We also analyze the 

' entropy and surface gravity for the horizon of the rotating dark matter solution. 

G\ : 

' Keywords: Dark matter, dark energy, dark fluid, energy conditions, surface 

^ . gravity. 

^ ' 1. Introduction 

' In a recent paper [1], referred to as Paper I, an exact solution of Einstein's field 

equations has been derived, describing dark matter, whose energy-momentum tensor 
admits a non-perfect fluid having the negative pressure and the dark energy equation 
of state parameter w = —1/2. The most striking property of this solution is that 
the metric of the solution describes both the background space-time structure and 
the dynamical aspects of the gravitational field in the form of the energy-momentum 
tensor. That is, the mass m of the solution plays the role of both the curvature of 
the space-time (non- flat) as well as the source of the energy- momentum tensor with 
Tab 7^ (non-vacuum) measuring the energy density and the negative pressure. This 
gives the information that when the mass m of the dark matter solution is set to zero, 
the space-time will become the flat Minkowski space with vacuum structure Tab = 0. 
In the case of Schwarzschild solution, the mass plays only the role of curvature of 
the space-time geometry and cannot determine the energy-momentum tensor. Thus, 
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the Schwarzschild solution is non-flat vacuum space-time with Tab = 0. Here the 
advantage of the dark matter solution over the Schwarzschild one is that it is non-flat 
and non-vacuum space-time metric. 

Nowadays, various possible models have been suggested for the composition of 
dark matter in [2-11] and references there in. The observations of luminosity-redshift 
relation for the type la supernovae [12-14] suggest that the missing energy should 
possess negative pressure p and the equation of state parameter w = p/p [11]. It is 
emphasized the fact that the approach in the derivation of the dark matter solution 
[1] is merely based on the identification of the power n = 2 of the Wang-Wu mass 
function [15] without any extra assumption. The identification of the power n = 2 in 
the Wang-Wu mass function has first considered in [1] and not been seen discussed 
before in the scenario of exact solutions of Einstein's field equations. This fact can be 
seen from the mass function as a power series of the coordinate r as 

+ 00 

M(n,r)= (ln{u)r'', (1.1) 

n=— oo 

where Qniu) are arbitrary functions of u. Wang and Wu analyzed the mass function for 
generalization of non-rotating Vaidya solution by choosing the function q-niu) corre- 
sponding to the number n [15]. In [16] the mass functions has been utilized in rotating 
system and found the role of the number n in generating rotating embedded solutions 
of Einstein's field equations. The roles of the power n in the expansion series (1.1) for 
the known spherically axisymmetric solutions are cited as follows [15, 16]) 

(i) n = corresponds to the term containing mass of the vacuum Kerr family 
solutions such as Schwarzschild, Kerr. 

(ii) n = — 1 purveys the charged term of Kerr family such as Reissner-Nordstrom, 
Kerr- Newman . 

(iii) n = 1 furnishes the global monopole solution [15]. 

(iv) n = 2 provides the dark matter solution [1] . 

(v) n = 3 contributes the de Sitter cosmological models, rotating and non-rotating 
[16]. 

These values of n are used for studying stationary solutions (non-rotating and rotat- 
ing). The non-stationary (time dependent) Vadya-Bonnor black holes, non-rotating 
[15] and rotating [16] can be obtained when n = and n = — 1 as Vaidya-Bonnor 
metric describes a charged solution. Non-stationary rotating as well as non-rotating 
de Sitter cosmological solutions can also be obtained when n = 3 [17]). In fact the 
Wang-Wu power series expansion of the mass functions happens to be the most conve- 
nient approach to generate new (embedded or non-embedded) solutions of Einstein's 
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field equations, using the complex scalar fields developed in [16] for the mass func- 
tion M{u,r,6) in Newman-Penrosc (NP) formalism [18]. The above identifications of 
power n indicate that the case n = 2 corresponds to the dark matter solution of Paper 
I, whose the line element is reproduced here for use as 

ds^ = {I - 2mr)du^ + 2dudr -r'^{de^ + siii^e dcfP') (1.2) 

where the constant m is considered to be the mass of the dark matter and is non-zero 
for the existence of the dark matter distribution. When u = constant, the surface is 
the future directed null cone. This line-element will not be asymptotically flat when 
r — > oo, and has a coordinate singularity at r = (2m) ~^ describing a Lorentzian 
horizon. The line element (1.2) is certainly different from Schwarzschild solution with 
guu = 1 — 2M/r having singularity at r = 2M, M being the Schwarzschild mass. The 
energy-momentum tensor for the above solution takes the following form 

Tab ={P + P){UaUb - VaVb) - PQab, (1-3) 

where Ua is a time-like vector {uau"" = 1) and is a space-like {vav"" = —1). gab is 
the metric tensor. The energy density p and the pressure p are found as 

4 2 

The parameter of the dark energy equation of state is the ratio of the pressure to the 
energy density and takes the constant value as 

It is noted that the energy-momentum tensor (1.3) having negative pressure (1.4) does 
not describe a perfect fluid. Hence, we refer it to as dark fluid with the equation of 
state (1.5). The 4-velocity vector fields Ua associated with (1.3) is expanding 6 7^ 0, 
accelerating ita 7^ as well as shearing Gab 7^ with zero- twist. This means that 
the observer of the dark matter does not follow the path of the time-like geodesic 
as Ua-bu'^ 7^ 0. The non-geodesic condition is a striking feature of the solution to be 
interpreted as an accelerated expansion of space-time geometry (1.2). The energy- 
momentum tensor for the dark fluid violates the strong energy condition p > 0, p + 
p > 0. The violation of strong energy condition is due to the negative pressure of 
the matter field content in the space-time geometry, and is not an assumption to 
obtain the solution (like other dark matter models mentioned in Sahni [2]). This 
violation indicates that the gravitational field of the solution is repulsive leading to 
the accelerated expansion of the space-time geometry [1]. The expansion of the dark 
matter solution with acceleration is consistent with the observational data [12-14]. 
Here it is to mention the fact that the energy-momentum tensor for the electromagnetic 
field in Ressner-Nordstrom and Kerr-Newman solutions does not violate the strong 
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energy condition p > 0, p + p > 0, since p takes positive value [16] causing attractive 
force of gravitational field, rather then repulsive as in the case of dark matter [1]. 

Prom the analysis of the dark matter solution in Paper I, we find that the dark 
energy density is only ascribed from the mass of the dark matter. It is the fact that 
without the mass of the dark matter one cannot determine the energy density and the 
negative pressure in the energy-momentum tensor distribution (1.3) in order to obtain 
the dark energy equation of state w = —1/2. This means that the negative pressure 
and the energy density of the dark fluid associated with the energy-momentum tensor 
are measured by the mass that produces the gravitational field in the space-time 
geometry of the dark matter. Hence, we may conclude that the solution (1.2) can 
explain the essential part of Mach's principle - "The matter distribution influences 
the space-time geometry" [19]. Here it is also emphasized that usually dark matter and 
dark energy are considered to be separate components of the Universe [2-11]. Arbey 
[20] has unified both components into a single "dark fluid" in order to explain the 
cosmological observation, based on a complex scalar field. In this regard the solution 
of Paper I indicates that dark matter and dark energy cannot be regarded as separate 
components, since the dark energy with the negative pressure is produced by the mass 
of dark matter, i.e. the same mass m plays the important key role for the existence of 
both. 

The space-time of the dark matter has an event horizon at r = (2m) Eventually, 
the area, entropy as well as surface gravity associated with the horizon have been 
analyzed. It is shown in [1] that the surface gravity is directly proportional to the 
mass of the dark matter. This indicates that the existence of the mass implies the 
existence of its surface gravity and the temperature of the dark matter on the horizon. 
The existence of dark matter horizon is in consistency with the cosmological horizon 
of de Sitter space with constant A [21], which is considered to be a common example of 
dark energy with the parameter w = —1 [3-11]. This parameter of equation of state is 
also satisfied for both the cosmological constant A as well as the cosmological function 
A{u) of the rotating and non-rotating de Sitter solutions [17]. It is quite interesting to 
specify that the approximate size of the dark matter mass is found, in the length scale 
r > 10^*^ suggested by Bousso [5], as m < (1/2) x 10"^'^, which is doubly bigger than 
the size of the cosmological constant |A| < 3 x 10~^^° with the horizon ta = a/3/A. 
This incredibly small size of the mass may be one reason why dark matter is not 
been able to observe, however its gravitational effect is being recognized in the form 
of dark energy in the dark halo of the Universe. Another possible reason may be that 
it behaves like a invisible black hole with event horizon. We believe that the non- 
rotating solution [1] may provide an example of an invisible dark matter which can 
produce dark energy with the equation of state parameter li; = — 1/2 in the accelerated 
expanding space-time geometry. However, the non-rotating dark matter solution of 
Paper I cannot explain the gravitational field of a rotating system. Therefore, this 
shows the requirement for developing a rotating dark matter solution of Einstein field 
equations in order to understand the complete nature of such invisible dark matter in 
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the Universe. 

Thus, the purpose of the paper is to generahze the rotating dark matter solution 
from the non-rotating one of Paper I, and then to investigate the nature of the ro- 
tating system. For instance, the rotating vacuum Kerr solution is an extension of 
the non-rotating Schwarzschild vacuum solution, and the Kerr-Newman solution is 
the charged rotating black hole which is an extension of the non-rotating Reissner- 
Nordstrom black hole. It is the fact that all astrophysical objects in the Universe 
are rotating on their own axes with very distinct physical properties. Hence, such a 
generalized rotating solution may provide more information about the gravitational 
field for better understanding the nature of dark matter in the Universe. 

The paper is organized as follows: In Section 2 the derivation of rotating dark 
matter solution of Einstein's equations of gravitational field has been shown. We find 
that the mass of the rotating dark matter solution proposed here describes the gravi- 
tational field of the space-time geometry and also determines the matter distribution 
with negative pressure, whose dark energy equation of state has the negative value 
w = —1/2 at the poles 9 = 7r/2 and 6 = 3it/2; otherwise it takes the functional 
form with minus sign as in (2.13) below. This is one of the curious properties of the 
rotating solution that the non-rotating one (1.2) cannot explain this situation. For 
the non-rotating dark matter a = it is always with the constant value it; = —1/2 as 
seen from the above. The variation of the value of w is consistent with those of the 
cosmological constant A [16] and function A{u) [17] of rotating de Sitter spaces. So 
we refer the solution obtained here to as rotating dark matter solution. The energy 
produced by this dark matter with the equation of state parameter is also regarded 
as rotating dark energy. However, the energy-momentum tensor of the dark matter 
distribution with negative pressure does not describe a perfect fluid, which can be seen 
in Section 2. This non-perfect fluid distribution is in accord with the remark of Islam 
[22] - "it is not necessarily true that the field is that of a star made of perfect fluid" . 
We also find that the dark matter solution has an event horizon. Consequently, we 
present the area, entropy and surface gravity for the horizon. The existence of the 
dark matter horizon discussed here is in accord with the cosmological horizon [21] of 
the rotating de Sitter spaces with constant A having the equation of state parameter 
10 = —1 at the poles. We conclude the paper in Section 3 with considerable remarks 
and evolution of the rotating solution. Thus, we summarize the results of the paper 
in the following theorem: 

Theorem 1 The rotating stationary dark matter solution is a non-asymptotic, non- 
vacuum, Petrov type D space-time, whose stress energy-momentum tensor possesses a 
dark fluid having a negative pressure and dark energy equation of state parameter with 
minus sign. 

Theorem 2 The 4-velocity vector of the rotating dark fluid having negative pressure 
is non-geodesic, expanding, accelerating, shearing as well as rotating. 
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Theorem 1 shows the physical interpretation of the solution that the space-time metric 
for the solution is not asymptotically flat when r ^ oo, and is Petrov type D with 
the non-zero component ip2 = —Cabcd^""'^^^'^'^'^ 7^ of Weyl tensor Cabcd- The 
energy-momentum tensor associated with the solution has a negative pressure with 
the dark energy equation of state parameter w{r,6) in (2.13). It is to emphasize that 
the Theorem 1 above is regarded as a generalization of the theorem of Paper I - 
"T/ie non-rotating stationary dark matter solution is a non- asymptotic, non-vacuum, 
conformally flat space-time, possessing a stress energy-momentum tensor of a dark 
fluid with negative pressure and dark energy equation of state parameter w = —1/2". 
Theorem 2 indicates the kinematic properties of the dark fluid content in the rotating 
axisymmetric space-time geometry of the dark matter. The presentation of this article 
is precisely based on mathematical calculations for obtaining a new rotating stationary 
(time independent) exact solutions of Einstein's field equations without any extra 
assumption. We do not attempt to use any observational data available at present. The 
Ncwman-Penrose (NP) formalism [18] in (—2) signature is adopted as mathematical 
tools in this paper. 

2. Rotating dark matter solution 

In this section we extend the non-rotating solution of Paper I to a rotating dark 
matter solution, as the rotating space-time metric has considerable potential appli- 
cation in astrophysical studies - all objects in astrophysics are rotating on their own 
axes. For instance, the non-rotating Schwarzschild vacuum black hole is generalized 
in the form of rotating Kerr vacuum one which has great application in the study of 
black holes with rotating configuration. Hence it is the main aim to derive a general- 
ized version of the non-rotating solution of [1] in the form of rotating stationary (time 
independent) axisymmetric dark matter solution for better understanding the nature 
of dark matter problem. For this purpose we choose the Wang-Wu function qn{u) in 
(1.1) as in Paper I 

, i m, when n = 2 

'^"^")= |o, when 2, ^^'^^ 

such that the mass function (1.1) takes the form 

+00 

M{u,r)= Yl qn{u)r'' = mr'', (2.2) 

n=— 00 

where m is constant and u = t — r is the retarded time coordinate. Then using the 
mass function (2.2) in general rotating metric presented in equation (6.4) of [16], we 
obtain a rotating metric, describing a stationary solution of dark matter in the null 
coordinates system {u, r, 6, (/)). The line element is found as 

ds^ = {1 - 2r^mR~^}du'^ + 2du dr + Aar^mR'^ sin^O du dcp 
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-2asm'^edrd(f> - R^dO"^ - {(r^ + a^f - Aahm'^ejR-hm'^e dcj)'^ , (2.3) 



where R"^ = + a^cos^O and A = — 2r^ m + a^. Here the constant a describes the 
rotational parameter. When a is set to zero, then the hne element (2.3) will reduce 
to non-rotating dark matter solution of [1]. This space-time metric is invariant to the 
simultaneous inversion of the retarded time u —u and the azimuthal angle (/) —(j) 
and also the metric coefficients be independent of u and (j). These are the requirement 
characters for a metric to be stationary and axisymmetric. The line element is non- 
asymptotic when r — oo. However, it has a singularity when A = having three roots 
given in (2.28) below. The line clement is certainly different from the Kerr metric with 
A* = — 2rM + o^, where M is the mass of the Kerr black hole having two roots for 
A* = 0. 

Then the covariant complex null tetrad vectors for the metric can be chosen as 
follows 

= 2^ -^a + - ^ « sin^^? (2.4) 

m„ = {ia sin e5l-R^5l- i{r^ + a^) sin 6 5^}. 

where R = r + ia cos 6. Here ia, Ua are real null vectors and nia is complex with the 
normalization conditions ia'n"' = 1 = —main"' and other inner products are vanished. 
These null tetrad vectors will be utilized to calculate Newman-Penrose quantities in- 
cluding spin coefficients, tetrad components of Ricci and Weyl tensors. From the tetrad 
components of Ricci tensor we calculate the energy- momentum tensor by virtue of Ein- 
stein's field equations. The Newman-Penrose spin coefficients [18], which are twelve 
complex quantities, representing the 24 Christoffel symbols involved in calculating the 
Ricci tensor are given below 

/{* = (j = A = e = z/ = 0, 
1 * A 



R' ^ 2RR'^' 
2ia — R cos 9^1 ^ 

(X = — , p = — r-^ cot &, 

2^2RRsm9' ' 2^2R 

iasinO iasinO 

The Ricci scalars of the metric are obtained as follows: 

4>ii = - Ji?a6(/V + m«m^) = ^^{r^ - ^a^os^O) 

A* ^ ^Rabg"^ = ^rm (2.6) 
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where Rab is the Ricci tensor. The complex tetrad component of the Weyl tensor Cabcd 
determining gravitational field of the space-time (2.3) is found as 

V'2 = -Cabcdl^m^rffvf- = -^^iar cos 9. (2.7) 

This expression of Weyl curvature scalar ip2 shows the coupling of the rotational pa- 
rameter a with the mass m of the dark matter. When the rotational parameter takes 
the limit a = 0, the Weyl scalar ■02 vanishes to lead to the confer mally flat space-time 

as in non-rotating dark matter solution of Paper I. From (2.7) we find that the rotating 
stationary dark matter (2.3) is type D in the Pctrov classification of space-time with 
a repeated principal null congruence £a (2.4), which is geodesic (k — 6 — 0), shear free 
((7 = 0), expanding {6 = ^^"a) as well as non-zero twist (a;^ = ^^[a;6]^"'^)) where 

^^^-li,'.^f = -^. (2.8) 

We find the energy-momentum tensor describing matter field for the above rotating 
stationary space-time (2.3) as follows: 

Tab = "^P^ianb) + 2pm(^arnb), (2.9) 
where the density and pressure are found as 

P = -17^2^{^ + Sa^cos^e}. (2.10) 



KE?B?' ^ KB?E? 

Because of the involvement of the rotational parameter a in a?cos^O, we find 

the difference in the values of p and p in (2.10) from those given in (1.4) of non-rotating 
solution. The energy-momentum tensor (2.9) can be expressed in the form of (1.3) 
using the orthonormal tetrad vectors {ua^Va^WaiZa} given in Paper I. Prom (2.10) it 
is found that the stress energy tensor (2.9) admits the weak energy condition {p > 0, 
p + p > 0) as well as dominant energy condition (p^ > 0, — > 0). However, the 
strong energy condition (p > 0, p + p > 0) is violated as the pressure p has negative 
value in (2.10). This violation may lead to a non-attractive (repulsive) gravitational 
force of a rotating dark matter solution (2.3) as pointed in [23]. 

The energy-momentum tensor of the solution (2.3) satisfies the energy conservation 
equations [17] formulated in Newman- Penrose formalism [18]. 

T"^^ = 0. (2.11) 

This establishes the fact that the space-time metric (2.3) describing dark matter is an 
exact rotating stationary solution of Einstein's field equations. The trace of energy 
momentum tensor Tab (2-9) is found as 

12r 

T = 2ip-p) = j^m. (2.12) 
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-^{r^ + 3a'^cos'^e) = 1. (2.14) 



Here it is observed that p — p > for rotating stationary of dark matter model. 
The energy-momentum tensor (2.9) does not describe a perfect fluid, i.e. for a non- 
rotating perfect fluid T^^^^ = {p + p)ua Uh — p gab with unit time-like vector Ua and 
trace T^p^^ = P — 3p, which is different from the one given in (2.12). 

The equation of state: The equation of the state for the dark energy is the ratio 
of the pressure to the energy density, and is found for the rotating solution (2.3) with 
(2.10) as 

w=- = -\\^ir'^ + ^a'^cos^e)\, (2.13) 
/9 2 I- J 

which shows the advantage of the rotating solution over the non-rotating one (1.2). 
Here we have seen the difference between the non-rotating (1.2) as well as the rotating 
dark matter solution (2.3) that the non-rotating one has the constant parameter w = 
— 1/2 of the energy equation of state; but that of the rotating solution takes the 
functional form (2.13), which may attend the constant value w = —1/2, if 

1 

This is true at the poles = 1^/2 and = 37r/2 with a 7^ 0, showing that the solution 
(2.3) describes the rotating dark matter producing dark energy with the equation of 
state parameter (2.13). This indicates the necessity for the analysis of a rotating dark 
matter providing more information than the non-rotating one [1]. It is noted that 
the equation (2.13) is always satisfied for the non-rotating dark matter solution of [1] 
with a = 0. However, the fact that w is not always equal to —1/2 for every value 
of 6 in (2.13), except 9 = 7r/2 and 37r/2, suggests the variation of equation of state 
w depending on the rotational configuration of the solution with a 7^ 0. It is the 
fact that we need not necessarily concentrate with non-rotating solution as in [1], and 
now we have a rotating solution which may provide the complete structure of dark 
matter in the Universe. The functional parameter w{r,9) (2.13) explains the main 
development of the rotating dark matter (2.3) producing dark energy, rather then the 
constant parameter w (1.5) for non-rotating solution. From (2.7), (2.9) and (2.13), we 
come to the conclusion of the proof of Theorem 1, stated in the introduction above. 

Kerr-Schild ansatz : Let us express the line-element (2.3) of the dark matter 
space-time in Kerr-Schild ansatz in order to show that it is certainly a solution of 
Einstein's field equations. For this purpose, we write the line element (2.3) as 

ds"^ = dt^ - dr^ - R^de^ - (r^ + a^) sin^ Q #^ - 2a sin^ Q dr d(p 

27T7 T"^ 

-^{du-asin^edcpy (2.15) 

R 

where dt = du+dr. This is of the Kerr-Schild ansatz on the Lorentzian flat background 

9T = 9ab-^^aib, (2.16) 
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with the null vector £a given in (2.4) and the Lorentzian fiat metric 

9af, = SlSi - 5151 - RHX - (r^ + a^) sin^ e5t5t - 2a sin^ ^<5^5^ (2.17) 
Now by the following transformation 

2 2 

dt = du+ ^ dr, d(j) = d(f) + ^dr, (2.18) 

we express the line-element (2.15) in Boyer-Lindquist coordinate system for future use 
as 

y\ sin^ 
ds^ = -^{dt - asrn^ e d(t)Y -^{{r"^ + a^)d(t) - adtY 

-^dr''-RV, (2.19) 

where '~' is being removed from d(p for simplicity. The line-element (2.19) can also 
be expressed in a very familiar form 

ds^ = -^{A - siii^ 9}df + ^^r^sm'^ 9 dtdcf) 

-^dr^ _ i?2rf02 _ 1 1(^2 + ^2)2 _ Aa^sm^e}sm^ed<P^. (2.20) 

where A = — 2mr'^ + a^. This form of the space-time metric is quite useful for 
working in the Boyer-Lindquist coordinate system {t, r, 9, (p). For example, we observe 
that the line element has a singularity when 

A = - 2mr^ + = 0. (2.21) 

This situation cannot be seen directly in (2.3). We can also find the stationary limit 
gtt = 0, which implies that 

- 2mr^ + cos 6 = 0. (2.22) 

Here we have seen that the singularity A = and the stationary limit gu = coincide 
at the poles 9 = 7r/2 and 9 = 3tt/2. One important point is to mention that the 
functional parameter w{r,9) (2.13) of the energy equation of state takes the constant 
value w = —1/2 at these poles. This situation cannot be observed in the case of 
non-rotating dark matter solution (1.2) of the Paper I. 

In order to express the rotating solution in Kerr-Schild anzatz based in Minkowski 
fiat background, let us use the following transformation of coordinates [24,25]. 

X = {r cos(f) + a sin^) sin^, y = {r sin^ — a cos<^) sm9, 

z = rcos9, t = u + r, (2.23) 

where r is defined in terms oi x,y, z as 

/ - (a;^ + y'^ + z'^ - a^)r^ - a^z^ = 0. (2.24) 
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The metric (2.3) can be expressed in the coordinate system {t, x, y, z) as 

ds^ = df - dx^ - dy^ - dz^ - M adx'^ f, (2.25) 

where the null vector (2.4) is being used to obtain 

iadx"' = dt — {r{xdx + ydy\ + a(xdy — ydx)} zdz. 

(H + a^) r 

The above transformed metric (2.25) is of the Kerr-Schild form 

gT=Vab + 2H{x,y,z)iaib, (2.26) 

where rjab is the Minkowski flat metric, g"^^ denotes the metric tensor for the rotating 
dark matter solution, and 

H{x,y,z) = - . , (2.27) 

is null with respect to both r]ab and y^™. The line element (2.25) is regular at every 
point except at r = 0, and z = 0. The two Kerr-Schild forms in different coordinate 
systems (2.16) and (2.26) establish the conclusion that the rotating stationary axisym- 
metric solution (2.3) is certainly a solution of Einstein's equations of gravitational field. 
By virtue of these Kerr-Schild forms, we may refer the rotating dark matter (2.3) to 
as dark matter solution is embedded into the flat background or simply dark matter in 
flat background. 

Surface gravity: Here let us show that the rotating dark matter metric has an 
event horizon at the point r = ri, which is a real root of equation (2.21) describing 
singularity of the solution. In fact the equation (2.21) has three roots ri, r2 and rs; 
ri is real and r2 and rs are complex conjugate of each other and are found as 

ri = ^{1 - h~ h-^}, (2.28) 



6m 

r2 = -^{l + ^{l + ^^/3) + ^{l-i^/S)h}, and = fs, 



where 



= {-1 - 54a^m^ + 6am^3(l + 27a'^m'^) } (2.29) 
These roots have the following relation 

(r — ri)(r — r2){i — ra) = {r^ — 2mr^ + a^j. 

2m 

We are interested only in the real root r = ri as the complex roots r2, r^ have 
unphysical meaning. When a = 0, we have h = —1 and the real root becomes ri = 
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(2m) ^, showing the consistency with that of the non-rotating solution (1.2). The 
area of the dark matter at event horizon r = ri is found as 



A = ^/mg^dOdcl) =47r{r? + a2}. (2.30) 

JO JO '■='■1 

According to Bekenstein-Hawking area-entropy formula S = A/4: [26], we find the 
entropy associated with r = ri as 

S = Tr{rl + a^}. (2.31) 

The gravity of the dark matter horizon is determined by the surface gravity which 
is defined as Kn"" = n^X/hn"" in [24]. Here the null vector n"' given in (2.4) above is 
parameterized by the coordinate u, such that d/du = n'^Va, and has the normalization 
condition iau'^ = 1 with the null vector The surface gravity k of the horizon 
associated with r = ri is found as 

« = -(7 + 7) = -4(l-3mr) . (2.32) 



r=ri 



where i?^ = + a^cos^^ and 7 is the spin coefficient given in (2.5). Here we do 
not intend to substitute the value of ri in (2.30-2.32) as these equations do not have 
simplified forms. Instead, we look at the following case of extreme dark matter. 

Extreme dark matter: Like extreme black holes [26], we shall discuss a case of 
extreme dark matter when 1 + 27a^m^ = in (2.29). In this case we have the quantity 
h = 1, which leads to the real root ri = — (6m)~^. Eventually, we find the area and 
entropy for the extreme dark matter as follows 

A = 7ra^, and S = Ja^, (2.33) 

where a is related with the mass m by virtue of the equation 1 -|- 27a^m? = 0. Then 
the surface gravity for this case is found as 

K = 27m(3-4cos^0)"\ (2.34) 

measuring directly from the mass m of the extreme dark matter, but can never be 
vanished for the existence of the dark matter. The functional parameter (2.13) of the 
energy equation of state becomes 

w = ~{l-4cos^e} (2.35) 

at ri = — (l/6m) of the extreme dark matter. This still retains the characteristic 
feature of the rotating dark matter having constant parameter w = —1/2 at the poles 
9 = 7r/2 and = 37r/2. 

Generalized Chandrasekhar's Theorem: Let us show the fact that the ro- 
tating dark matter solution (2.3), which is non- vacuum {T^b 7^ 0) Petrov type D 
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space-time, satisfies the condition of generalized Chandrasckhar's theorem [27]. Chan- 
drasekhar [25] has estabhshed a relation of spin coefficients p*, fi* , r, vr in the case 
of an affinely parameterized geodesic vector, generating an integral which is constant 
along the geodesic in a vacuum {Tab = 0) Petrov type D space-time 

P U T IT 

^ = ^ = - = -. (2.36) 

P* fl* IT T 

The original derivation of this relation is necessarily based on the vacuum Petrov type 
D space-time with Weyl scalars ^"2 7^ 0, ■00 = '01 = V's = '04 = and vanishing of 

the tetrad components of Ricci tensor 0oi = 4>02 = 0io = 'P20 = 4>i2 = 02i = 0oo = 
022 = 011 = A = 0. In [27] a complex tetrad component xi = ^fab{f'""n^ + fn°"mJ^) of 
Killing- Yano tensor fab has been introduced in the above relation (2.36) as 

C = C = - = - = -^> (2-37) 

p* p* TT T Xl 

where fab is the skew-symmetric tensor f^ab] satisfying the Killing- Yano equations 

fab;c + fac;b = 0. (2.38) 

The importance of KY tensor in General Relativity seems to lie on Carter's remark- 
able result [24] that the separation constant of Hamilton-Jacobi equation (for charged 
orbits) in the Kerr space-time gives a fourth constant. In fact, this constant arises 
from the scalar field Kabv"" which has vanishing divergent along a unit vector 
tangent to an orbit of charged particle. Here Kab = fmafb^ is a symmetric Killing 
tensor. The KY tensor is skew-symmetric of rank two and has six real components. 
These six components can be expressed by three complex tetrad components xoj Xii 
X2- Due to the vanishing of Wcyl scalars = 0i = '03 = "04 = for a Petrov type 
D space-time ip2 7^ 0, the tetrad components X2 of Killing- Yano tensor become 
zero. Then the remaining tetrad component xi can be solved from the Killing- Yano 
equations (2.38) transcribed in NP formalism in [27], and is found as 

Xl = iC{r - ia cos9) (2.39) 

where C is a real constant. Using this in (2.37) we obtain the relation as 



p1 = = L = 1 = = ^ 
p* p* Tt f Xl 



(2.40) 



where R = r + ia cos9 and the spin coefficients p* , p*, r, vr are given in (2.5) for the 
dark matter solution. Then the Killing- Yano tensor for the solution is found 

fab = 4a cos e n[„4] -I- Air m[amb] , (2.41) 

and accordingly, the Killing tensor Kab = fma fj^ becomes 

Kab = -8{a^ cos^ ^ + r^rn(^amb)}. (2.42) 
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We observe that the non- vacuum rotating dark matter solution is one of the examples 
which satisfy the generalized Chandrasekhar's relation (2.37). The non- vacuum ro- 
tating stationary solutions, like Kerr-Newman, Kerr-Newman-de Sitter, negative mass 
naked singularity embedded into de Sitter space [28] and the dark matter solution dis- 
cussed here are all Petrov type D space-times and admit the relation (2.37). Hence, 
we conclude, without any lost of generality, that the generalized relation (2.37) with 
Killing- Yano scalar xi holds true for every rotating stationary Petrov type D solution, 
vacuum or non-vacuum, admitting the KY (2.41) and Killing (2.42) tensors. 

Properties of time-like vector: Here we shall discuss the nature of the 4- velocity 
vector u"", which is time- like, appeared in the energy-momentum tensor for the dark 
fluid, if expressed in terms of orthonarmal tetrad vector fields {ua,Va,Wa, Za} as in 
(1.3), for the rotating dark matter solution (2.3) with (2.10). In fact the 4-velocity 
vector fields explain the kinemetic properties of a fluid distribution whether the 
fluid flow is expanding (0 = u^.^j 7^ 0), accelerating {iia = Ua-^u^ 7^ 0), shearing 
c"ab 7^ or rotating {wah 7^ 0). For example, the dark matter of Paper I is expanding, 
shearing, accelerating but non-rotating. The key role of the rotational parameter a of 
the stationary solution (2.3) can be seen in the following quantities. The expansion 
scalar B and acceleration vector iia = Ua-^^v!' are found as 

© ^ <a = ^^(l + 3mr2), (2.43) 
2r 2 A -1 a^sin^cos^ - 

where Va = -^{^a — ng) and A = — 2mr^ + . Then the shear Uah and vorticity 
Wah tensors are respectively obtained as follows 

msin6' _ . iasin^ . ^\ - ,^ . .\ 

- ^^^2 (^'^ + cos 0)v^amb) - ^j^j^2 ~ ^)v{amb) , (2.44) 

Wab = U[a-b] - U[aUb] 

mcos 9 , , , , , 

+71^^ ^ ^"'^"""'^' ^ ^ ^ 

which are orthogonal to {i.e., (TabU^ = and WabU^ = 0). However, the mass m 
directly determines the expansion G as well as the acceleration Ua of the solution as 
seen in (2.43). The existence of the vorticity tensor is due to the presence of the 
rotational parameter a in (2.45), and conflrms that the solution (2.3) is a rotating 
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stationary solution with accelerated expansion (2.43). When a sets to zero, the space- 
time will become the non-rotational solution (1.2) with vanishing the vorticity tensor. 
The existence of the acceleration vector shows that the observer of the rotating dark 
matter solution discussed here follows the non-geodesic path of the time-like vector 
{ua:bU^ 7^ 0). This provides the key feature for the accelerated expansion of the rotating 
dark matter. Prom Equations (2.43)-(2.45) imply the proof of the Theorem 2 stated 
in the introduction above. 

3. Conclution 

In this paper we derive an exact stationary (time independent) axisymmetric so- 
lution of Einstein's field equations describing rotating dark matter whose energy- 
momentum tensor possesses a non-perfect fluid having the negative pressure and the 
non-constant parameter w{r^ 6) of equation of state with minus sign. However, it takes 
a constant value w{r,0) = —1/2 at the poles 9 = 7r/2 and 9 = 37r/2. The most ex- 
iting property of this solution is that the metric of the solution describes both the 
background space-time structure and the dynamical aspects of the gravitational field 
in the form of the energy-momentum tensor. That is, the mass m of the rotating 
solution plays the role for both the curvature of the space-time (non-flat) as well as 
the source of the energy-momentum tensor with T^b ^ (non-vacuum) measuring the 
energy density and the negative pressure (2.10). In the case of rotating Kerr solution, 
the mass plays only the role of curvature of the space-time and cannot determine the 
energy-momentum tensor. Hence, the Kerr solution is non-flat vacuum space-time 
with Tab = 0. Here the advantage of the solution (2.3) over the Kerr is that it is a 
non-flat and non-vacuum rotating stationary space-time metric . 

The expansion, acceleration, shear and vorticity of the 4-velocity vector fields Ua 
suggest that the observer of the rotating dark matter follows the non-geodesic path as 
Ua-fiU^ 7^ 0- We find that the energy-momentum tensor for the rotating solution violates 
the strong energy condition p > 0, p + p > due to the negative pressure (2.10) of the 
matter field content in the space-time. This violation indicates that the gravitational 
field is repulsive leading to an accelerated expansion of the rotating universe. This 
expansion of the dark matter solution with acceleration (2.43) is consistent with the 
observational data [12-14]. 

It is emphasized the fact that the approach of the derivation of rotating dark matter 
solution here is merely based on the identification of the power n = 2 of the Wang-Wu 
mass function without any extra assumption. This identification of the power n = 2 
in the mass function (1.1) has considered first in Paper I as mentioned above. It is to 
mention here that Kerr metric can be obtained by choosing n = in the Wang-Wu 
mass function (1.1). In the analysis of the rotating dark matter solution, we find that 
the dark energy density is only attributed from the mass of the dark matter. It is the 
fact that without the mass of the dark matter we cannot measure the energy density 
and the negative pressure of the energy-momentum tensor distribution in order to 
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obtain the dark energy equation of state w with minus sign. This means that the 
negative pressure and the energy density of the dark fluid associated with the energy- 
momentum tensor arc measured by the mass that produces the gravitational field in 
the space-time geometry of the rotating dark matter solution. The most important 
property of the rotating stationary dark matter solution (2.3) is that the parameter of 
the dark energy equation of state becomes a function w{r,9) (2.13) with minus sign. 
The functional parameter w{r^6) provides constant value —1/2 at the poles 9 = tt/2 
and 9 = 37r/2, retaining the important character of a dark energy even in the case 
of rotating system. These poles are important points for the rotating dark matter 
that the event horizon and the stationary limit are coincide at these poles, where the 
functional parameter of the equation of state takes the constant value w{r, 9) = —1/2. 
This situation could not be observed in the case of the non-rotating dark matter 
solution of Paper I. The rotating solution (2.3) suggests that dark matter and dark 
energy cannot be regarded as separate components, since the dark energy with the 
negative pressure is produced by the mass of the dark matter, i.e. the same mass m 
plays the important role for the existence of both. 

The rotating space-time metric appears singular when A = 0, which has three 
roots, one is real and other two are complex. The real root r\ corresponds to an event 
horizon. The event horizon of the dark matter is very similar to the one surrounding 
a black hole. Accordingly, we find area, entropy as well as surface gravity for the 
horizon. The existence of dark matter horizon is consistent with the cosmological 
horizon of de Sitter space with constant A [21], which is considered to be a common 
candidate of dark energy with the parameter w = —1 [3-11]. This parameter w = —1 
of the equation of state is also true for both the cosmological constant A as well as the 
cosmological function A(7i) of the rotating and non-rotating de Sitter solutions [17]. 

From the analysis of the dark matter solutions we find the profound changes that 
(i) the non-rotating dark matter solution (1.2) is non- asymptotic, non-vacuum, con- 
formally flat metric; (ii) the rotating stationary solution (2.3) is a non- asymptotic, 
non-vacuum, Petrov type D in the classification of space-times, whose one of the re- 
peated null direction la is geodesic, shear-free, expanding [9 = \l"'-a — \{P* + P*)^] 
as well as non-zero twist [Cj"^ = ^£[a;b]^"''^ = -{{p* + P*)^] such that Cabc[d'^k]'^''£'^ = 0, 
where p* is the spin coefficient given in (2.5). It is noted that to the best of the authors 
knowledge, the solution describing rotating stationary dark matter presented here is 
not been seen discussed before. The embedded version of the rotating solution into 
Kerr black hole as dark matter in Kerr black hole background may be seen elsewhere 
in another paper. We believe that this rotating solution may provide an example of 
an invisible rotating dark matter, which can produce dark energy having the equation 
of state parameter with minus sign in an accelerated expanding space-time geometry. 
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